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We present a theoretical study of the conductance and thermopower of single-molecule junctions 
based on Ceo and Ceo-terminated molecules. We first analyze the transport properties of gold-Ceo- 
gold junctions and show that these junctions can be highly conductive (with conductances above 
O.lCo, where Go = 2e^ /h is the quantum of conductance). Moreover, we find that the thermopower 
in these junctions is negative due to the fact that the lowest unoccupied molecular orbital (LUMO) 
dominates the charge transport, and its magnitude can reach several tens of /xV/K, depending on 
the contact geometry. On the other hand, we study the suitability of Ceo as an anchoring group 
in single-molecule junctions. For this purpose, we analyze the transport through several dumbbell 
derivatives using Ceo as anchors, and we compare the results with those obtained with thiol and 
amine groups. Our results show that the conductance of Ceo-terminated molecules is rather sensitive 
to the binding geometry. Moreover, the conductance of the molecules is typically reduced by the 
presence of the Ceo anchors, which in turn makes the junctions more sensitive to the functionalization 
of the molecular core with appropriate side groups. 
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I. INTRODUCTION 

The fullerene Ceo is attracting a lot of attention in the 
field of molecular electronics.^ One reason is that the de- 
localization of the frontier orbitals of Ceo suggests that 
this molecule can be a good candidate to build highly 
conductive single- molecule junctions, a goal that remains 
elusive and that it has only being achieved with short 
molecules. 1^^ On the other hand, it has been recently 
suggested that Ceo used as an anchoring group to bind 
molecules to the electrodes can improve the reproducibil- 
ity of the conductance measurements in single-molecule 
junctions,^ which is a crucial issue in molecular electron- 
ics. The goal of this work is to further analyze these two 
different questions from a theoretical point of view. 

The first experiment on individual Ceo molecules was 
reported by Joachim et al^ Here, a scanning tunneling 
microscope (STM), with a tip made of tungsten, was 
used at room temperature to study the conductance of a 
Ceo molecule on an Au(llO) surface. It was shown that 
in the contact regime this hetero junction has a conduc- 
tance of 2.35 X 10~^Go, which is clearly lower than in 
the case of metallic atomic-size contacts. Since then, dif- 
ferent groups have investigated experimentally the trans- 
port properties of Ceo molecular junctions, mainly with 
gold electrodes, and they have reported very different re- 
sults. Thus for instance. Park et a/.l^ performed measure- 
ments in Au-Ceo-Au junctions using the electromigration 
technique and depositing the Ceo molecules from a liquid 
solution. In this case, the conductance at low bias was 
found to be largely suppressed and the current-voltage 
characteristics were dominated by the Coulomb blockade 
phenomenon. A related experiment, but this time with 
a microfabricated break junction, showed a much higher 



low-bias conductance (of the order of O.lGo), which was 
attributed to the appearance of Kondo physics. ^^ A more 
systematic study of the low-bias conductance was carried 
out by Bohler et al}^ln this case, the authors conducted 
low-temperature (10 K) break junction experiments in 
which the molecules were evaporated in situ. From the 
analysis of a conductance histogram, it was concluded 
that Au-Ceo-Au junctions exhibit a preferred conduc- 
tance value close O.lGo- On the contrary, recent STM 
break junction experiments at room temperature have 
shown a very large spread of conductances with a certain 
preference for values around 5 x 10~^Go-^^ 

The main evidence that Ceo junctions can have a 
rather high conductance has been provided by a se- 
ries of controlled STM experiments with other electrode 
materials (different from gold) performed by Berndt 
and coworkers and which have been nicely backed up 
theoretically.^^^Hi^ Thus for instance, Neel et al}^ re- 
ported a controlled STM study in ultrahigh vacuum 
(UHV) in which Ceo molecules deposited onto copper sur- 
faces exhibited conductance values of the order of 0.25Go 
in the contact regime. It is also worth mentioning that 
Kiguchi^^ reported break junction experiments at room 
temperature in UHV conditions in which a single Ceo 
molecule between Pt electrodes was shown to exhibit a 
conductance as high as O.TGo- 

Although the electronic and transport properties of 
Ceo metal- molecule- metal junctions have been addressed 
theoretically by numerous groups, ^^"^^ ^^ ^^ studies of 
transport in Ceo junctions with gold electrodes, the most 
commonly used metal in molecular electronics, are sur- 
prisingly rather scarce. ^^ ^^ Moreover, those references 
explore only ideal contact geometries which may have 
little to do with those realized in the experiments. For 



this reason, we present here an ab initio study of the 
charge transport in Au-Ceo-Au junctions paying special 
attention to the role of the contact geometry. Our analy- 
sis, based on the combination of density function the- 
ory (DFT) and nonequilibrium Green's function tech- 
niques, shows that conductances above 0.1 Go are possible 
in realistic contact geometries. Moreover, motivated by 
very recent experiments,^^ we have investigated the ther- 
mopower of these junctions and found that this quantity 
is negative, as reported in Ref. 12, which is simply due to 
the fact that the low-bias transport is dominated by the 
Ceo LUMO. Furthermore, we have found that the ther- 
mopower varies significantly with the contact geometry, 
and its magnitude can reach several tens of /iV/K, which 
is higher than in previously investigated molecules. ^^^^ 

The second topic that we want to address in this work 
is the role of Ceo as an anchoring group. Some of the main 
challenges in the field of molecular electronics are related 
to the fabrication of single- molecule junctions with very 
well-defined transport properties and the ability to tune 
those properties at will. A strategy that is being pursued 
to achieve these goals is the use of suitable anchoring 
groups to bind the molecules to the metallic electrodes. 
The thiol (— SH) group is the most commonly used an- 
choring group, specially when the electrodes are made 
of gold, because of their high covalent bond strength.!^ 
However, the thiol group has been shown to lead to a 
large variety of binding geometries,^^^^^ which implies a 
large spread in the observed conductance values. Many 
different alternatives to the thiol group have been ex- 
plored in recent years. For instance, Venkataraman and 
coworkers^^ introduced the amine group (— NH2) as an 
interesting possibility to obtain better defined values in 
the conductance histograms, which was attributed to a 
higher selectivity of amine-gold binding. Similar con- 
clusions have been drawn in a recent analysis of nitrile- 
terminated (— C=N) biphenyls.'^Sl The list of anchoring 
groups explored in molecular junctions increases steadily 
and the search for the "most" convenient group, leading 
to highly reproducible transport properties, has becom e 
one of the central issues in molecular electronics. ^^^ ^^ 

In this context, Martin et alJ put forward the inter- 
esting idea of using Ceo as a new anchoring group. The 
idea is that Ceo offers a large contact area which, to- 
gether with the high molecular symmetry, may reduce 
the spread of the conductance values. Moreover, this 
fullerene is known to strongly hybridize with metallic 
surfaces,^ and, as explained above, it has been shown in 
different STM and break junction experiments that it can 
sustain a rather high conductance. Indeed, in Ref. 7 the 
electrical characteristics of l,4-bis(fullero[c]pyrrolidin-l- 
yl)benzene (BDC60) (with Ceo anchor groups) were stud- 
ied using gold microfabricated break junctions, and it 
was found that the conductance histograms exhibited 
more pronounced peaks than those obtained with 1,4- 
benzenediamine and 1,4-benzenedithiol. More recently, 
Leary et al.^^ have shown that the use of Ceo as anchor- 
ing group facilitates enormously the characterization of 



single-molecule junctions in STM experiments under am- 
bient conditions and it allows to establish unambiguously 
the conductance of the molecule under study. 

These experimental results are promising, but it re- 
mains to be explored whether the use of Ceo as a ter- 
minal group still allows, for instance, for the possibil- 
ity to chemically tune the conductance by an appro- 
priate functionalization of the molecular core, a s it h as 
been demonstrated with other anchoring groups IMIMSI] 
In other words, the main role of an anchoring group must 
be to provide the chemical link to the electrodes without 
modifying the essential properties of the molecular back- 
bone. In this sense, it remains to be shown whether or 
not Ceo is too invasive to be used as an anchoring group. 

Besides the results for the pure Ceo, we also present a 
study of the transport properties of molecular junctions 
based on Ceo-terminated molecules: BDC60 and several 
derivatives. Our DFT-based analysis aims at addressing 
two main questions: (i) Does Ceo reduce the spread in 
conductance values found with other anchoring groups? 
and (ii) Is Ceo too invasive to be used as a suitable an- 
choring group? Our results suggest that the conductance 
and thermopower of Ceo-terminated molecules are still 
quite sensitive to the binding geometry and we expect a 
large spread of values in typical STM and break junction 
experiments. On the other hand, our results indicate that 
Ceo may reduce the electron communication between the 
molecular core and the metallic electrodes, leading to a 
reduction of the conductance. However, this reduction of 
the effective metal-molecule coupling and the fact that 
the frontier orbitals lie relatively close to the Fermi en- 
ergy lead to a notable increase in the sensitivity of the 
junctions to the functionalization of the molecular back- 
bone, as compared with thiol or amine groups. 

The rest of the paper is organized as follows. In the 
next section we briefiy describe the methodology em- 
ployed to compute the transport properties of single- 



molecule junctions. Then, in Sec. |III| we present a de- 
tailed analysis of the conductance and thermopower of 



Au-Ceo-Au junctions. Sec. IV is devoted to the analysis 
of junctions with BDC60 molecules modified by the in- 
clusion of several side groups. The results for the conduc- 
tance and thermopower of these junctions are compared 
with those obtained using thiol and amine as anchoring 
groups. Finally, we summarize the main conclusions of 
this work in Sec. [Vl 



II. METHODOLOGY 

Our main goal is to describe the transport proper- 
ties of single- molecule junctions based on Ceo molecules 
and Ceo-terminated compounds. For this purpose, we 
employed the DFT-based transport method described 
in detail in Ref. [^, which is built upon the quantum- 
chemistry code TURBOMOLE G.l.E^In this method, the 
first step is the description of the electronic structure of 
the molecular junctions within DFT. In all our calcu- 



lations we used the BP86 functionaP^ and the def-SVP 
basis set.ESlln order to construct the junction geometries, 
we first relaxed the molecules in the gas phase. Then, the 
molecular junctions were constructed by placing the re- 
laxed molecules between two finite clusters of 20 (or 19) 
gold atoms and performing a new geometry optimization. 
In this optimization, the molecule and the 4 (or 3) out- 
ermost gold atoms on each side were relaxed, while the 
other gold atoms were kept frozen. Subsequently, the size 
of the gold clusters was extended to around 63 atoms on 
each side in order to describe the metal-molecule charge 
transfer and the energy level alignment correctly. 

The final step in our method is to transform the in- 
formation about the electronic structure of the junctions 
obtained within DFT into the different transport prop- 
erties. This is done using nonequilibrium Green's func- 
tion techniques, as described in detail in Ref. 021 In 
the coherent transport regime, and following the spirit 
of the Landauer approach, the low-temperature linear 
conductance is given by G = Gqt{Ey) = Go X^i '^i(^F), 
where Go = 2e^ /h is the quantum of conductance, t{Ey) 
is the junction transmission at the Fermi energy, £^f, 
and {Ti{E)} are the transmission coefficients, i.e., the 
energy-dependent eigenvalues of the transmission matrix. 
The second transport property of interest in this work is 
the thermopower, which within the coherent transport 
regime is given by 
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with Kn{T) = j dE{E - iiYT{E)[-dEf{E,T)], where 
ji is the electrochemical potential and f{E^T) = [1 + 
ex.-p[{E—/jj)/kBT]]~^. We shall compute the thermopower 
at room temperature (T = 300 K), and in all the exam- 
ples discussed here one can still use the low-temperature 
expansion of Eq. (IT]) , which is given by 



TT^klTr'jEF) 
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Here, the prime denotes a derivative with respect to en- 
ergy. Thus, the thermopower measures the logarithmic 
first derivative of the transmission function at £^ = E^. 
The sign of this quantity carries information about the 
location of the Fe rmi energy within the gap of a molecu- 
lar junction.'^^l^ 



III. CONDUCTANCE AND THERMOPOWER 
OF GOLD-Ceo-GOLD JUNCTIONS 

This section is devoted to the analysis of the trans- 
port properties of Au-Ceo-Au junctions, which will also 
serve us as a reference for the study of the Ceo-terminated 
molecules. Let us start our analysis by recalling the elec- 
tronic structure of Ceo in the gas phase. Within our 
DFT approach, and in agreement with Ref. [24l we find 
that the highest occupied molecular orbital (HOMO) is 
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FIG. 1. (Color online) (a) Two ideal Au-Ceo-Au junctions 
with top (left) and hollow (right) binding geometries, (b) The 
transmission as a function of energy for the two geometries 
in (a), (c-d) In these two panels the solid line corresponds 
to the total transmission, while the others correspond to the 
contribution of the individual transmission coefficients as a 
function of energy. 



five- fold degenerate, while the LUMO is three- fold de- 
generate (at —5.90 and —4.26 eV, respectively). These 
energies have to be compared with the Fermi energy of 
gold, which in our calculations is —5.0 eV. In order to elu- 
cidate how the electronic transport takes place through 
a Ceo molecule coupled to gold electrodes, we first con- 
sider two ideal geometries in which the molecule is bound 
to the electrodes in a top and in a hollow position, see 
Fig. Ilia). These geometries have been constructed as 
follows. We first relaxed the molecule on top of a single 
cluster, then we added a second cluster symmetrically 
at the other end, and finally we relaxed again the whole 
junction, as described in the previous section. In the top 
position, we find that the apex gold atom binds to two 
carbon atoms of a 6:6 bond, each C-Au distance being 
around 2.45 A. This geometry is consistent with that re- 
ported in Ref. [21] for various Ceo-gold nanocontacts. In 
the hollow position, similar to that explored in Ref. [23l 
the three- Au-atom terrace is facing one carbon atom and 
C-Au distances are in the range 2.3-2.4 A. 

In Fig. fllb-d) we show both the total transmission 
and its channel decomposition as a function of energy 



for these two ideal contact geometries. The first thing to 
notice is that the conductance, determined by the trans- 
mission at the Fermi energy, is very high as compared 
with other organic molecules of similar length, 0.55Go 
and 1.85Go for the top and hollow position, respectively. 
In both cases the low-bias conductance is dominated by 
the LUMO of the molecule, as it has been found in STM 
experiments of Ceo on Au surfaces (see e.g. Ref . l48|) . For 
the top position, we find that the transmission at the 
Fermi energy is largely dominated by a single channel 
which originates from one of the LUMOs of the molecule, 
which is split from the other two and is shifted to lower 
energies due to its better coupling to the electrodes. For 
the hollow-type geometry, we find that two conduction 
channels give a significant contribution to the low-bias 
conductance. These channels originate from two of the 
LUMOs, which in this case are more strongly coupled to 
the electrodes than in the top geometry due to the higher 
number of C atoms in direct contact with the electrode 
atoms. This is simply the reason for the higher conduc- 
tance of this geometry, which agrees with the findings of 
Ref. ITS] where it was shown (both experimentally and 
theoretically) that the conductance of a Ceo junction in- 
creases with the number of atoms in contact with the 
molecule. Let us also mention that conductances above 
IGo have also been reported in theoretical studies of Ceo- 
junctions with Al,^^ ^^ Au,^^ and Cu^^ electrodes when 
the leads are similar to ideal surfaces, i.e., with a high 
Au-Ceo coordination. 

The calculations discussed above suggest that Au-Ceo- 
Au junctions can have a conductance comparable to 
metallic atomic-size contacts. However, the ideal geome- 
tries considered so far should provide a rough estimate 
for the expected conductance values since it is unlikely 
to realize experimentally contacts with such a high de- 
gree of symmetry. Thus, a more direct comparison with 
the experiments requires a detailed analysis of the junc- 
tion formation and of the evolution of the conductance 
during the stretching of the contacts. This is precisely 
what we have done, as we now proceed to explain. In or- 
der to simulate the junction formation, we started with 
a geometry in which the molecule is positioned laterally 
with respect to the gold-gold axis and we used gold clus- 
ters terminated with a single Au atom. Fig. [2ja). Then, 
the gold electrodes were separated step-wise (in steps of 
~ 1 A) and the junction geometry was relaxed in ev- 
ery step. This protocol was repeated until the junction 
was broken and the molecule lost contact with the elec- 
trodes. To characterize the junction during the stretch- 
ing process, we computed different quantities such as the 
binding energy of the junction, the Mulliken charges in 
the Ceo molecule, the linear conductance, and the ther- 
mopower at room temperature. The results are shown in 
panels (b) to (e) of Fig. ^ Notice that in panel (d) we 
have also included the conductance of a Au-Au junction 
with the same distance separation to estimate how much 
current is fiowing directly from gold to gold (bypassing 
the molecule) in the different stages of the elongation 
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FIG. 2. (Color online) (a) Some representative geometries of 
the stretching simulation of a Au-Ceo- Au junction. They cor- 
respond to Au-Au distances (distance between the Au tips) 
of 4.2 A, 10.2 A, and 17.5 A. The other panels show the 
following quantities during the stretching process: (b) bind- 
ing energy of the junction, (c) charge on the Ceo molecule, 
(d) conductance of the Au-Ceo- Au junction (circles) and, for 
comparison, conductance of a Au-Au junction with the same 
distance separation (squares), and (e) thermopower at room 
temperature. 



process. 

In our simulation, after the first steps, the molecule ro- 
tates and then it places itself in the middle of the junction 
adopting a geometry in which the top gold atom is bound 
to a single C atom, see central geometry of Fig.[2Va). This 
structure, which is the most stable one with a binding en- 
ergy close to 1 eV, differs from the very symmetric top 
geometry in Fig.[TJa). In this geometry the Au-Ceo inter- 
action is maximized by the proximity of the side surface 
of the Au cluster. This is consistent with a related anal- 
ysis reported in Ref. 20", Then, after further stretching, 
the center of the molecule is aligned with the junction 
axis and the molecule remains there until the rupture 
of the contact. In the Au-Au distance range 8-12 A, 
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FIG. 3. (Color online) Transmission as a function of energy 
for the Au-Ceo-Au junctions shown in Fig.l3|a). The different 
curves correspond to different elongation stages, as indicated 
in the legend {Az refers to the Au-Au distance). 



where the binding energy is maximum in magnitude, the 
molecule is negatively charged and the conductance ex- 
hibits a "plateau" with values between 0.07 and 0.2Go, 
which is consistent with the experiments of Ref. TT, At 
an Au-Au distance of ~12 A, the contact breaks, as sug- 
gested by the evolution of the binding energy, and the 
conductance starts to decrease exponentially, while there 
is a tiny positive charge in the molecule. 

The thermopower results shown in Fig. |2|e) deserve 
special attention in view of the recent experimental re- 
sults reported in Ref. [121 In that work, thermopower 
measurements of fullerene- metal junctions were per- 
formed at room temperature with a STM break junc- 
tion technique. In particular, for Au-Ceo-Au junctions a 
preferred value of —14.5 ± 1.2 /iV/K was found, the mi- 
nus sign suggesting that the LUMO dominates the con- 
duction. As we show in Fig. [3| where we display the 
transmission curves of Au-Ceo-Au junctions at different 
elongation stages in the simulation, the low-bias trans- 
port is dominated by the LUMO of the molecule at any 
distance. As a consequence, the thermopower is negative 
at any stage of the elongation process, see Fig.|2|e), and 
in particular, its value for the most stable geometry is 
approximately —35 /iV/K, which is a factor two larger 
than in the experiment. Let us also mention that for 
the ideal junctions of Fig.jlja) we obtain a thermopower 
of —19.37 /iV/K for the hollow geometry and a value of 
—91.62 /iV/K for the top one. This big difference be- 
tween these two geometries is due to the fact that the 
transport takes place in an almost on-resonant situation. 

To conclude this section, it is worth commenting that 
in our simulations we have not taken into account the 
van der Waals interactions. In this sense, the binding 
distances might not be exact. Dispersion forces play a 
role in the interaction between Ceo and gold."^ However, 
the Au-Ceo binding is known to be mainly covalent with 
some ionic character j^^l^^l^^* This is corroborated by the 
fact that, during our simulated elongation, the molecule 



is pulled in between the electrodes, due to the chemical 
interaction. 



IV. Ceo AS AN ANCHORING GROUP 

We now analyze the role of Ceo as an anchoring group 
in molecular junctions. We have seen in the previous 
section that this molecule can sustain a rather high 
conductance, which suggests that Ceo can, in princi- 
ple, provide a very efficient electronic communication, 
when used as an anchoring group. To explore this 
idea, let us first study the transport through a 1,4- 
bis(fullero[c]pyrrolidinyl) benzene (BDC60) molecule, see 
Fig. [4J in which a phenyl ring is connected to two 
fullerenes on two opposite sides via a pyrrolidine group 
(in a so-called "dumbbell" fashion). We have chosen this 
molecule for several reasons. First, it has b een in vest i- 
gated both experimentall}^ and theoretically ]^ ^^^^ -^ which 
allows us to establish a comparison with our results. Sec- 
ond, the transport through the central moiety (a phenyl 
ring) can be analyzed with other anchoring groups, which 
is necessary to determine the quality of Ceo as a terminal 
group. 

Our DFT calculations of the electronic structure of the 
isolated BDC60 molecule show that its HOMO appears 
at -4.7 eV and that it is localized on the central part 
of the molecule, while the two- fold degenerate LUMO is 
localized on the Ceo's, as displayed in Fig. [4J These two 
LUMOs are the lowest in energy of a series of six levels 
(ranging from -4.18 to -3.87 eV) which originate from the 
interaction between the three LUMO orbitals of each Ceo, 
as explained in detail in Ref. 52, In agreement with this 
reference, we find that the nitrogen atoms are displaced 
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FIG. 4. (Color online) (a) HOMO and (b-c) two-fold degen- 
erate LUMO of the BDC60 in the gas phase. 
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FIG. 5. (Color online) Total transmission and individual 
transmission coefficients as a function of energy for the two 
Au-BDC60-Au junction geometries shown in the upper part 
(top and hollow positions). 



from the phenyl ring plane. 

Turning now to the analysis of the Au-BDC60-Au junc- 
tions, we have again studied the conductance of two dif- 
ferent types of geometries, a hollow and a top binding 
geometry, see the upper part of Fig. |5J In the hollow 
position, the three top Au atoms are bound to four C 
atoms, while in the top position, the apex gold atom is 
bound to two C atoms of a 6:6 bond. In Fig. [5] we show 
the results for the total transmission and the channel de- 
composition as a function of energy for both geometries. 
In both cases, the transmission close to the Fermi energy 
is dominated by a single channel and the resonance just 
below Ef originates from the HOMO of the molecule. 
In spite of the fact that the HOMO is pinned very close 
to the Fermi energy in both cases, the conductance is 
equal to 9.0 x 10~^Go foi" the hollow position, while it is 
2.5 X 10~^Go for the top geometry. These values should 
be compared with the preferential value of 3 x 10~^Go 
reported in the experiments of Ref. |7l although a large 
spread of conductance values was also found there. We 
attribute the low values of the conductance (as compared 



TABLE I. HOMO and LUMO energies of the molecules (in 
gas phase) based on the phenyl unit and different anchoring 
and side groups. 



to the Ceo junctions) and the difference between the two 
geometries to the weak effective coupling of the phenyl 
ring to the Ceo molecules, which is quite apparent in 
the small width of the transmission resonances. In other 
words, the phenyl- Ceo effective coupling is the actual bot- 
tleneck in these junctions and its weakness makes the 
conductance very sensitive to the exact level alignment 
and to the metal-molecule coupling. 

Let us compare our results for the BDC60 molecule 
with other theoretical results published recently. First, 
we find that the current is mainly carried by the HOMO 
of the BDC60, while in Ref. |52] it was found that the 
transport is dominated by the LUMO. Let us stress that 
we have confirmed the level alignment described above 
by test calculations with even larger gold clusters (116 
atoms). The discrepancy between these results may be 
due to differences in the electrodes' shape (in Ref. 52 the 
electrodes were modeled as ideal surfaces) and to the pe- 
riodic boundary conditions applied in their modeLl^ Sec- 
ond, in Ref. 52 it was claimed that the conductance of the 
Au-BDC60-Au junctions is not very sensitive to the bind- 
ing geometry, while we find a large difference between the 
top and hollow geometries. We attribute this discrepancy 
to the fact that in that reference no binding with un- 
dercoordinated Au atoms was considered. On the other 
hand, in Ref. 51 it was stated that the weak coupling and 
the insufficient conjugation throughout the three parts of 
this molecule is detrimental for the electronic transmis- 
sion. Although our results can not be directly compared 
with those of Ref. [51] (supporting LUMO transport), we 
agree in observing that the peaks corresponding to the 
frontier orbitals (HOMO and LUMO) appear narrower 
and lower than in the transmission curve of Ceo- 

After studying BDC60, we want to address the is- 
sue of whether or not the conductance of this dumb- 
bell molecule can be chemically tuned by functionaliz- 
ing the phenyl unit, as it is known to be possible with 
other anchoring groups. To this aim, we have investi- 
gated three different substituents: CH3, F, and CI, and 
we shall compare the results with those obtained employ- 
ing two other widely used anchoring groups, namely thiol 
(— SH) and amine (— NH2). The CH3 group is known 
to be electron-donating, while F and CI are electron- 
withdrawing groups, and indeed for the molecules with 
thiol and amine anchoring groups the HOMO and LUMO 
are pushed upward in energy when the molecule is func- 
tionalized with CH3, while they are pulled downward 
with F and CI, see Table [l[ However, in the case of the 




FIG. 6. (Color online) Geometries of the studied molecular 
junctions with a phenyl ring functionalized with chlorine and 
attached to the gold electrodes via (a) thiol, (b) amine, and 
(c) Ceo groups. 



dumbbell molecules this trend is not reproduced. The 
LUMO is not affected by the presence of the side groups, 
as expected since it resides in the Ceos, and the HOMO 
is shifted to lower energies also in the presence of the 
CH3 group, see Table [ij Moreover, in this case the func- 
tionalization causes a distortion of the central part, due 
to the interaction between the pyrrolydine and the sub- 
stituents. This distortion, which does not occur in the 
case of SH and NH2 because of their lower steric hin- 
drance, is responsible for the unusual behavior of the 
CH3 side group. 

As an illustration, in Fig|6]we show the molecular junc- 
tions for the three different anchoring groups and with 
CI functionalization. We choose a top binding geometry 
for NH2 (let us remind that the amine group only binds 
to undercoordinated Au sites) and a hollow one for SH 
and Ceo- In Fig[7|we show the transmission curves for 
all molecules with different anchoring groups, as well as a 
comparison of the conductance values in the lowest panel. 
The first thing to notice is that for all anchoring groups, 
the relative energy positions of the frontier orbitals upon 
side functionalization reproduce the trends observed in 
the gas phase. For thiol- and amine-terminated molecules 
the current flows through the HOMO, consistent with 
what was found for tolane molecules with the same an- 
choring groups. f^ Concerning the conductance values, 
the functionalization has no dramatic effect in the cases 
of the thiol and amine group, see lower panel of Fig [7[ 
The shift in the position of the frontier orbitals is more 
apparent in the values of the thermopower, which are 
shown in Table [ill for all the molecules. In this table 
one can see the conflrmation of the naive expectation 
that says that when the transport is dominated by the 
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FIG. 7. (Color online) From top to bottom: Transmis- 
sion curves for molecules anchored via thiol, amine, and Ceo 
groups, and conductance values for all the junctions. 



HOMO, the thermopower increases when the HOMO is 
shifted to higher energies, and it decreases when this or- 
bital is pushed to lower energies. Moreover, it is worth 
stressing that the thermopower for the benzenedithiol 
molecules has been measured by Baheti et alW-" and our 
results are in good quantitative agreement. 

As one can see in Fig [71 the effect of the side groups is 
much more pronounced in the case of the Ceo-terminated 
molecules. In particular, the functionalization in this 
case lowers considerably the conductance. The main 
reason for that is that the transmission resonance that 





Thermopower (/iV/K) 


CH3 


H 


F 


CI 


s 
NH2 

Ceo 


8.36 
8.96 
18.90 


7.31 
6.98 
96.74 


5.18 
6.29 
13.0 


4.03 

5.64 

-12.95 



TABLE II. Thermopower of the different molecular junctions 
analyzed in Fig. ^ with different side groups and anchoring 
groups. 



dominates the transport, and which is associated to the 
HOMO of the molecule, is much narrower for this an- 
choring group and therefore, it is much more sensitive to 
the shift induced by the side group. Notice also that in 
the case of the CI group, the conductance is indeed dom- 
inated by the LUMO of the molecule due to the strong 
energy shift of the HOMO in this case. This fact is re- 
flected in a change of sign in the thermopower (see Ta- 
ble [ll|, which is something that does not occur for the 
other two anchoring groups. 



V. CONCLUSIONS 

In summary, we have presented a DFT-based analysis 
of the conductance and thermopower of individual Ceo 
and Ceo-terminated molecules with gold electrodes. We 
have shown, in agreement with several experiments, that 
Au-Ceo-Au junctions can have a rather high conductance 
above O.lGo for realistic geometries. Moreover, we have 
found that the transport through Ceo takes place through 
its LUMO, which leads to a negative thermopower in 
agreement with recent measurements. The fact that the 
LUMO lies relatively close to the Fermi energy, which 
means in practice that the energy derivative of the trans- 
mission at the Fermi energy is rather large, leads to a 
rather high thermopower in comparison with other or- 
ganic molecules. 

On the other hand, to investigate the use of Ceo as 
an anchoring group, we have first studied the trans- 
port through Au-BDC60-Au junctions and found that 
the conductance is rather sensitive to the binding geom- 
etry. Furthermore, we have found that the conductance 
is decreased, as compared with the Ceo junctions, due 



to the poor electronic communication between the Ceo's 
and the molecular core (phenyl unit). Then, in order to 
study whether Ceo as a terminal group is too invasive, 
we have analyzed several BDC60 derivatives which differ 
in the presence of a side group in the phenyl unit (CI, F, 
and CH3), and we have compared the results with those 
obtained using thiol and amine anchoring groups. Our re- 
sults indicate that the BDC60-based junctions are much 
more sensitive to the functionalization, i.e. the changes 
in the conductance and in the thermopower induced by 
the side groups are much more significant in the case of 
the molecules with Ceo as anchoring group. 

So in short, our study supports the idea that Ceo is 
a good conductor and it suggests that it can be used 
as a convenient anchoring group to study typical effects 
related to the chemical modification of the molecules: 
role of side groups, degree of conjugation, length de- 
pendence, etc. However, Ceo does not seem to resolve 
the usual problem related to the spread of conductance 
values. Moreover, in dumbbell molecules like BDC60, 
beside the substituent-related shifting effect, also config- 
urational changes due to steric repulsions can play an 
important role. 
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